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Abstract 

We present a first-principle path integral Monte-Carlo (PIMC) study of the binding energy of excitons, 
trions (positively and negatively charged excitons) and biexcitons bound to single-island interface defects in 
quasi-two-dimensional GaAs/Ala;Gai_a;As quantum wells. We discuss in detail the dependence of the bind- 
ing energy on the size of the well width fluctuations and on the quantum-well width. The numerical results 
for the well width dependence of the exciton, trions and biexciton binding energy are in good quantitative 
agreement with the available experimental data. 
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I. INTRODUCTION 



Excitonic atoms and molecules in quantum confined semiconductors have been intensively 
investigated in the last decade. These systems show nontrivial Coulomb correlation effects leading 
to interesting optical and transport characteristics not seen in bulk materials^ A strong increase of 
the binding energy of the excitonic complexes was found experimentally y, 
with decreasing quantum well (QW) width and increasing magnetic field | 

In the literature there has been an active discussion about the influence of localization potentials 

n 

on the bindin g en erg y of excitons and excitonic complexes HUH . Most of the theoretical calcula- 
tions 13 Qll^ Q [3. [3] show a substantially weaker dependence of the binding energies 



nji 




on the QW width than those experimentally observed lll|,|2l|4l|5l|6l|7|,l8l|9lH- In particular, it 
was found experimentally that the binding energy of the charged excitonic complex is for narrow 
quantum wells much larger than the one estimated theoretically. The difference being typically a 
factor of two for narrow QWs. An explanation for this could be the trapping of the excitons (trions, 
biexcitons) by ionized donors in the barriers lisl or by some kind of interface defects produced by 
the mixture of well and barrier materials during the QW growth process, i.e. by QW width fluc- 



tuations or fluctuations in the a 

inRefs. iy,y,y,y,y,n 



loy composition of the barrier, which were not taken into account 
|. Such effects can induce an additional weak lateral confinement 
which leads to the confinement of the particles in all three dimensions like in the case of a quantum 
dot potential. The low-temperature photoluminescence of such structures originates from the ra- 
diative recombination of the exciton states localized at such nonuniformities of the heterostructure 
potential. In this situation broadening and splittings of both exciton and trion peaks develops in 
the PL spectra ligj . Such lateral confinement becomes more important in narrow quantum wells. 
Theoretical calculations of exciton, trion and biexciton states in such structures is a fairly complex 
problem because of the need to take simultaneous account of the Coulomb interaction and the 
three-dimensional heterostructure potential, which is no longer translationally invariant. 

The standard theoretical approach to calculate binding energies is to solve the correspond- 
ing many-particle Schrodinger equation by means of an appropriate basis expansion. This works 
efficiently in simple geometries but is not easily applicable to our problem with well-width fluctu- 
ations. Recently, a different approach was developed which is based on solving the Bloch equation 
for the many-particle density matrix [24]. It was demonstrated in Ref. ll24ll that this problem can 
be efficiently solved using path integral Monte Carlo (PIMC) methods without any restrictions 
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on the geometry of the confinement potential. No quantum well width fluctuation effects were 
considered in Ref. 112411 . 

The aim of the present paper is to understand and explain recent experimental data on the 
binding energy of the ground state excitons, trions, biexcitons in QW's by including localization 
effects. We consider localization as a consequence of the local modulation of the thickness of the 
quantum well of 1 — 2 monolayers which corresponds to the experimental findings of Ref. Iigj . 
In agreement with experimental results, we find that such QW width fluctuations can increase the 
trion binding energy in GaAs-based quantum wells by up to 100 percent as compared to ideal 
QW's without interface roughness. 

We also found that for lateral localization diameters exceeding D 150 A the binding energy 
of the negative trion can become larger than that of the positive trion, in contrast to the case of ideal 
QW's where the positively charged excitons are slightly more strongly coupled. The reason is that 
the localization confinement has a different influence on the lateral wave functions of electrons 
and holes. Thus the trion composition (i.e. versus X ) becomes crucial to the value of the 
binding energy when the localization diameter changes. 

Our numerical method is an extension of the Path Integral approach of Ref. j24|. The method 
does not involve expansions in terms of basis functions, no symmetry assumptions are made (in 
this sense it can be considered as first-principle), and the error can be managed in a controllable 
way |25i]. 

The paper is organized as follows. In Sec.|ll|we present and discuss the Hamiltonian for the 
exciton, biexciton and charged excitons in a quantum well with interface defects. We also discuss 
the approximations used in the present calculations. In Sec. |ni| we introduce the basic ideas of 
the numerical method, i.e. Path Integral Monte Carlo (PIMC), used to obtain the ground state of 
the excitonic complexes. In Sec.|lV|we compare the correlation, localization and binding energy 
of the localized exciton (X), biexciton (X2) and charged excitons {X^) ground state with the 
ones of the non localized, i.e. free exciton complexes in a quantum well. Further, we study the 
dependence of the X, X^ and X2 ground state properties on the defect width and height. In 
Sec.|V|we compare our calculations with the available experimental data and, finally, present our 
conclusions in Sec.lVIl 
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II. THE THEORETICAL MODEL 



We consider a single GaAs quantum well grown between two Al^.Gai_^As barriers. The effec- 
tive mass framework is used to describe the semiconductor material and the QW structure. Using 
the isotropic approximation the Hamiltonian for N^. electrons and holes reads: 

i=l 

where rrij and are the mass and charge of the i-th particle, e is the dielectric constant, which we 
assume equal for the well and for the barrier, VJ= (^h) is the confinement potential associated with the 
presence of the QW, is the lateral (localization) confinement which is due to the fluctuations 
of the QW width. We take the quantum well growth direction as the z-direction. 

For a GaAs/Al^^Gai-^ As quantum well, we consider the following heights of the square-well 
potential: K = 0.57 x (1.155s + 0.37a;2) el^ for electrons and \4 = 0.43 x (1.1552; + O.STa;^) eV 
for holes. In our calculations we use an Al concentration of x = 0.3. Furthermore, the following 
material parameters are used: e = 12.58, mg = 0.067 ttiq, my, = 0.34 mg, where uiq is the mass 
of the free electron. The units for energy and distance are H* = 2R* = e^/ (e a^) = 11.58 meV, 
aB = h^e/{me e^) = 99.7 A, respectively. We have also considered the case of an anisotropic hole 
mass according to using for the in-plane hole mass a smaller value of m[J = 0.112mo, and 
in the quantum well growth direction ml = 0.377mo. Comparing the binding energies calculated 
with the isotropic and anisotropic approximations gives important insight about the relevance of 
band structure details for the properties of excitonic complexes in quantum wells. 

The actual shape of the interface defects is not known and depends on the sample growth condi- 
tions. To limit the number of parameters, we simulate the interface defects through a cylindrically 
symmetric potential with a lateral radius R and height V^°^. The potential height is determined by 
the zero-point energy and was obtained as the difference between the lowest energy levels of the 
electron (hole) in two QW's with the widths L and L + S, where 5 = na with n = 1, 2 and a is 
the thickness of a single monolayer. Because of the difference in mass between the electron and 
the hole, the height of the localization potential will also be different (see Fig.[T]). For large well 
widths (L) the localization potential is given by {f?-n'^ lmi)b jL^ . Notice that in Fig. [T]for small 
L the electron localization has a local maximum which is due to the increased penetration of the 
electron wave function into the barrier material. For the hole this occurs at much smaller L due to 
its larger mass. 
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FIG. 1: Dependence of the height of the locahzation potential V^°l^y Eq. for electrons (open dots) and 
holes (full dots) on the well width L for a well width fluctuation of one monolayer. Inset: the lowest energy 
level in the potential, V^i^f^y vs. the well width. 

In GaAs, 1 (2) monolayer(s) correspond to a well width fluctuation of 5 = 2.8 (5.6) A. These 
parameters ensure that the exciton (trion, biexciton) wave function in the growth direction z is 
practically not affected by the defect. It is very instructive to see from Fig. 1 that, for narrow QW's, 
this lateral localization potential reaches about 15 meV which is comparable to the exciton binding 
energy and is several times larger than the trion binding energy. This behavior is in qualitative and 
quantitative agreement with the monolayer splitting measured experimentally ligj . In the inset of 
Fig. [U we plot the value of the lowest energy level in a QW as a function of the QW width. The 
main figure can be obtained directly from the results of the inset through V''°^ = Eq {L+5)~Eq (L) . 

We proceed further with the assumption that the QW confinement is sufficiently strong and that 
the Coulomb interaction among the particles in the z-direction will not modify the wave functions 
in the ^-direction, and consequently we may use in the ^-direction the noninteracting electron and 
hole wave functions. In this adiabatic approximation we neglect the influence of in-plane electron- 
hole correlations on their motion perpendicular to the QW plane. This assumption is valid due to 
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the strong quantization in square wells of widths L < aB, giving rise to the condition on the 
energies: AE^^^^-^ ^ Ec, E^ , E^^,E^'^, where AElfj^-^ is the level spacing in the quantum well, 
and Ec, E^ , E^^ ,E^^ are the correlation and binding energy of exciton, biexciton and trions, 
respectively. 

Our approach to compute the binding energies starts from the A^-particle (A^ = 2, 3, 4) density 
matrix of the excitonic complex of interest (exciton, trion, biexciton) which is obtained from a 
solution of the corresponding Bloch equation, see Ref. 11241] . In the adiabatic approximation, the 
full A^-particle density matrix factorizes into 

where (R^^) = {r^i, r^s, • • • , v^n^; Yhi, r/,2, • • • , ^hNh) is a ?>D (2D) vector of all particle 
coordinates, ^e(h) is the z coordinate of all electrons (holes), p{Ze, {3) and p{Zh., (3) are the density 
matrices of free electrons and holes confined in the z direction by the square well, and /5 = 1/ fc^T 
is the inverse temperature. We underline that the density matrix p(R^^,/3) contains all in-plane 
electron-hole correlations and fully includes the effect of the localization potential. It obeys the 
two-dimensional A^-particle Bloch equation which is obtained by averaging the three-dimensional 
Bloch equation over z and using Eq. dU) and the ansatz in Eq. (O 



i=l 



Here, we have introduced an effective 2D in-plane interaction potential Vjf': 

K^yiP)= I dz^dz,, Yl i"^^' . piz,,p)piz,,p) 



X 



dZ^ dZh piZe, 13) p{Zh, 13) 



(4) 



and the total localization potential 



yioc,.y\E^(L + 5)-E,{L), if./WTf)<D/2- 

[o, if^{x^ + y^)>D/2, 

where Eq{L) is the lowest energy level in a QW of widths L (see inset of Fig. 1). 



ni. NUMERICAL SIMULATION APPROACH. 

We numerically solve the Bloch equation ^ using the path integral representation of the den- 
sity matrix. Using the operator identity e~^^ = (e"""^)^^, the density matrix at inverse tem- 
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perature P can be expressed in terms of M density matrices, each taken at a M times higher 
temperature M ksT, or as a path integral with M steps of size r = l/{MkBT) llisll . 



p(R, R; /5) = JdRi...J dKM-i J2 

(R|e-^^|Ri) . . . (RAf-i|e-^^|PR), (6) 



where R represents a set of coordinates of particles in 2 dimensions; P is the -particle ex- 
change operator, (— 1)"^^ denotes the sign of the permutation for Fermi particles (electrons and 
holes); p(R, R';r) = (R|e^^^|R') is the coordinate representation of the A^-particle density 
matrix at the new inverse temperature r. For the A^-particle high-temperature density matrix, 
p(R, R'; r), we use the pair approximation which is valid for r < 1/(3 H*), 

p(R, R'; r) - n p''' (r^' <^ ^) H ^"/^ ? t''' ^ + 

i fj-^^ P^^^ (i-i, K; t) pill (rfe, r'^; r) 

where i, j are particle indices and pl^^ (p'^O is the one(two)-particle density matrix. The one- 
particle density matrix, pl^^, is the known free-particle kinetic energy density matrix. The pair 
density matrix pl^^ was obtained from a direct numerical solution of the two-particle Bloch equa- 
tion for which we used the matrix squaring technique IitH isI. 

As one can see from Eq. Q, the needed diagonal matrix elements of the low-temperature den- 
sity operator are expressed in terms of all diagonal and off-diagonal matrix elements of the corre- 
sponding high-temperature density operator which can be effectively computed using path integral 



Monte Carlo simulations, see Ref. ll24ll and references therein. Obviously, for these simulations to 
be efficient, it is crucial that the off-diagonal density matrix, pl^^, can be quickly evaluated for any 
given initial (r^, r^) and final (r/, r^') radius vectors of the particle positions. For this reason, 
before doing the PIMC simulations, we calculated in advance tables of the pair density matrices 
(DM) for each type of interaction in our system. In our electron-hole system in a QW with the 
localization potential, we needed to calculate: i) three tables of pair density matrices correspond- 
ing to electron-electron, hole-hole and electron-hole interactions given by the two-particle Bloch 
equation with the smoothened effective 2D Coulomb potential, see Eq. Q, and ii) two tables with 
the density matrix of a single particle (electron or hole) in a 2D cylinder of finite height (for par- 
ticles localized at the interface defect). The contributions of all these interactions (correlations) 
can be treated as additive, once the used high-temperature pair density matrices correspond to suf- 
ficiently high temperature (such that commutators of pairs of energy contributions are negligibly 
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small). Finally, using the pair DM tables, we are able to calculate the many -body density matrix, 
Eq. ©, for any set of initial R and final R' positions of all particles. We substitute this expression 
into Eq. © and perfoj;m the high dimensional integration using a multilevel (bisection) Metropolis 
algorithm (see, e.g. 

In the present calculations we used tables of the pair density matrices at a temperature three 
times the effective electron-hole Hartree, i.e. 1 /r = 3H* = 403 K. By choosing in Eq. ^ the 
number of factors equal to M = 270, the full density matrix, p(R, R; /3) and all thermodynamic 
quantities can be accurately evaluated at a temperature T = 1.49 K. All results shown below 
correspond to this temperature value. 

Before considering in detail the effect of quantum well width fluctuations on the binding ener- 
gies of excitonic complexes, we recall the main results obtained for ideal QW's with finite width 
L 1I24I1 . Quasi-two-dimensional systems like GaAs QW's have been extensively investigated in the 
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last years, both experimentally yj. 
These studies revealed that, due to the confinement, the 2D excitonic states have binding energies 
which are several times larger than the binding energies in the bulk materials. This effect is mainly 
due to the confinement of the carrier wave functions along the structure growth direction, which 
leads to a two-dimensional character of excitons and, consequently, to a change in the in-plane 
interaction potential between the carriers. In the framework of the adiabatic approximation these 
changes can be easily seen in the effective in-plane potential , Eq. Q, which depends on the 
quantum well width L as a result of the integration in Eq. Q over the free electron and hole den- 
sity matrices which reflect the probability distributions of an electron and hole in the square well 
potential of given width L. In Fig.|2l we present the electron and hole densities in the square well 
with a width varying between 10 A and 160 A. The results in the left panel of Fig. |2l confirm that, 
due to the smaller mass, the electron is less localized than the hole and, for L < 20 A, most of the 
electron density resides in the barrier material. 

The effective in-plane potential is shown in the right panel of Fig.|2l Notice that it depends 
on L in a non-monotonic way reaching a maximum (absolute) value around L ^ 40 A. Such 
an increase with decreasing L (for not too small values of L) is also found experimentally and 
theoretically and is due to an increase in the interparticle correlations and results in the main 
contribution to the increase of the binding energies in ideal quantum wells at intermediate QW 
widths. 
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FIG. 2: Left: electron Pe{z) and hole ph{z) density matrix in the QW (dotted lines indicate QW walls). 
Right: QW width dependence of the effective electron-electron (ee) and electron-hole (eh) potentials, see 
Eq. ©. 

IV. BINDING ENERGIES 

In this section we investigate the combined influence of the finite QW width and of the interface 
defects (defect width and height) on the ground state of the exciton and excitonic complexes. In 
particular, we analyze the modification of the binding energies and of the average interparticle 
distances in the ground state of excitons (X), positive and negative trions (X^) and biexcitons 

(^2). 
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A. Binding energies and average size of excitonic complexes 



For an ideal QW, i.e. without interface defects, we define the binding energy of the exciton, 
charged exciton and biexciton as: 



where (/i) is the energy of a single electron (hole) in the given quantum well with a free particle 
mean kinetic (thermal) energy /c^T, and E(A) is the total energy of the excitonic complex A. If 
an interface defect is present and a localization potential is included in our calculations, then the 
above definitions must be modified. All energies must be replaced by the corresponding energies 
of particles localized in the defect potential. The corresponding generalized expressions will be 
given in Sec. IIVBI 

Using a finite temperature approach such as PIMC, one calculates states in thermal equilibrium. 
Moreover when the temperature is not sufficiently low and comparable with the depth of the trap- 
ping potential, the equilibrium state reached in a sufficiently long simulation will correspond to 
non-localized states rather than localized ones. To correctly obtain the total and binding energies 
of localized excitonic complexes, the results were computed not by averaging over all states, but 
by restricting the average to the states localized in the trapping potential. 

We now discuss the results for the binding energies and average interparticle distance in the 
ground state of various excitonic complexes as a function of the depth and width of the interface 
defects. In Fig.|3ta) we plot the binding energies versus the diameter of the trapping potential, D, 
for the case of a 1 monolayer (1 ML) surface defect. The corresponding relative gain in the binding 
energies due to the interface defect is shown in Fig. |3tb). As an example, we took a QW width 
of L = 60 A and a IML QW width fluctuation which corresponds to the following heights of the 
electron and hole localization potentials, \ V^°'^\ = 3.43 meV and \ Vl°''\ = 1.28 meV, respectively. 
From Fig. 121 one can notice that for all excitonic complexes the binding energy is always larger 
when a defect is present than in the ideal case. In particular, it increases with the diameter D of the 
trapping potential up to some maximum after which it slowly decreases. However, for very large 
D the system approaches very slowly the ideal QW result, but 1 ML wider than the original one. 

Notice that the position of the maximum is different for the different excitonic complexes. This 
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(8) 
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FIG. 3: (a): Binding energies of various excitonic complexes vs. the diameter of the 1 ML quantum well 
width fluctuation for a QW width of L = 60 A and temperature, T = 1.5 K; (b): The same as (a) but now 
for the relative increase of the binding energy. 

is readily explained by the different lateral size of different bound states which is determined by the 
lateral extension of the electron and hole wave functions in the trap and by their relative distance. 
The electrons are more sensitive to the defect because the trapping potential has a larger effect on 
their localization (the holes are substantially localized even in the absence of the defect). Further, 
we observe that the lateral confinement has a very different effect on the magnitude of the exciton, 
trion and biexciton binding energies, see Fig.|3tb). In particular, the exciton binding energy is only 
relatively weakly affected by the localization, i.e. a very small peak in the relative binding energy 
gain of less than 20% occurs. In contrast, the binding energy of the negatively charged exciton 
increases by more than 100%, from 1.4 meV to 3 meV for the localization potential of diameter 
D ^ 300 A. 

In Fig. m we show the average in-plane interparticle distance, pij, versus the diameter of the 
localization trap for the exciton, trions and biexciton. For a 2D system this results, after using the 
adiabatic approximation, in the following expression: 

rij g{rij) dvij / J gij-ij) dvij, (9) 

where g{rij) is the pair distribution function of the particles i and j. 

By comparing the electron-hole distances in various complexes, see Fig. H^b), it can be seen 
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that the electron-hole distance in the exciton, i.e. the size of the exciton, is about 1.2 — lA times 
smaller than the electron-hole separation in the charged excitons and about 3 times smaller than 
the average electron-electron (hole-hole) distances, Fig.Hta). This explains our previous finding, 
see Fig. |3b), that the exciton state is much less influenced by the lateral confinement than the 
X^. In the exciton, where the electron and hole are coupled much stronger, the interparticle 
distance changes only slightly with the diameter D and the effect on their binding energy is weak. 
Notice also that the peak in the gain of the binding energy quite closely follows the minima of 
the electron-hole interparticle distances. This result agrees well with the experimental findings 
(discussed below) that in the case of localized particles the binding energy of the exceeds the 
one of the X+. 

Further, it is interesting to note that the biexciton appears to be less extended than the trions, 
thus explaining the fact that trions have a lower binding energy than the biexciton, see Fig. Ha). 
This is consistent with the experimental observations, which will be discussed in detail in Sec.fVl 
At the same time, the biexciton is more affected by the interface defect than the positive trion, see 
Fig-Etb). This suggests that the number of electrons in the excitonic complexes plays a much more 
important role in the interaction with the interface defect than the number of holes. In fact, both, 
the X^ and X2, which contain two electrons, are more influenced by the localization than the X 
and X+. The reason is that the localization potential has a stronger impact on the confinement of 
electrons than on holes, as noted above. 

Next we compare the negative and the positive trion. For small localization islands, D < 150 A, 
the average distances between electrons in X~ and holes in X+ are very similar and, consequently, 
the relative gain in the binding energies of the two trions is close as well, see Fig.|3tb). In contrast, 
for wide localization islands, i.e. D > 200 A, the behavior of the two trions differs significantly, 
e.g. the binding energies and the interparticle distances between the respective pair of equally 
charged particles deviate from each other, see Fig.|3land Fig.|4| The reason is a different influence 
of the QW potential on the electron and hole wave functions in the X~ and X+ states. In a 
quantum well the electrons are substantially extended into the barrier material (i.e. in the z- 
direction), whereas the holes are much more confined, see Fig.|2ta). This leads to a significantly 
larger overlap of the two electrons in the X^ state compared to the two holes in the X+ and to a 
weaker in-plane effective electron-electron interaction potential, Eq. Q. As a result, in Fig.Ua) 
at D = we can see that the electron-electron distances, p^e, are slightly smaller than those of 
the holes, phh- Though the average distances in the X~ and X+ states are very similar, if one 
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FIG. 4: Average distance between the constituents of the different excitonic complexes as a function of D 
for (a) equally charged and (b) oppositely charged particles. 

compares fluctuations of the average distances, 5p = a/ (p^) — (p)^, they will be much larger for 
the electrons in the negative trion . This can be seen directly from the behavior of the pair 
distribution functions shown in Figs.|7ta,b). When the trions are localized by the additional lateral 
confinement these fluctuations are quenched and the distance between particles of the same charge 
is decreased, e.g. for < D < 300 A. This must have a stronger effect on the electrons than on 
the holes. The stronger repulsive interaction between the holes in the X+ state reduces to a larger 
extent the gain in the binding energy as compared to the state, where the electron interaction 
is weaker and particles can be brought to smaller distances. This accounts for the increase of the 
X~ binding energy up to, and even beyond, the X^ binding energy for large D. 



B. Origin of the enhanced binding energies 

It turns out that the differences in the trion binding energies are not caused solely by their 
different spatial extensions, see Fig. HI In fact, the difference between the average interparticle 
distances in the two trions is not sufficiently large to account for the gain of the binding energy 
for the X and the X+, see Fig. Elb). The explanation must be found in the fact that the gain 
in the binding energy, as a function of the size of the localization island, does not come only 
from the changes in the Coulomb interaction related to the interparticle distances, but also from 



13 



changes (and differences) of the electron and hole localization energies, and the kinetic energy of 
the particles. 

In the presence of the localization potential each electron and hole acquires an additional (single 
particle) potential energy - the localization energy -Efj^^ . This single particle energy can be directly 
computed by averaging the localization potential over the radial electron (hole) distribution, 

Kc'^(^) = / rfre(.)^7^(re(,))\/A°'^(r,(,))/ J dr g\ e^h)) ■ (10) 

Similarly, each bound excitonic complex is affected by the localization potential where each elec- 
tron and hole contributes additionally to the total localization energy 

eux) = ei^{x) + eUx), (11) 

where E'ij^^\x) for the exciton is computed in an analogous way as was done for the single 
particles in Eq. dlOb . but with the appropriate radial electron (hole) distribution inside the localized 
exciton. Further, Eq. (flOb can be straightforwardly generalized to the trions and biexciton cases. 
Obviously, the localization energy modifies the total energy of all particles, 

E{X) = EiX, D = 0)^ EiX, D)- E{X^) ^ E(X±, D)- EiX^) -> ^(Xa, D), (12) 

where for all bound states (X, X^^ X2) the total energy can now be written as 

E{D) = E(0) + 6Eco.x{D) + 6E^,{D) + Ei,,{D). (13) 

Here 5Ecou\{D) and 5E^n{E>) denote, respectively, the change of the Coulomb and kinetic energy 
due to the presence of the defect of diameter D. From Eqs. © and ( fT?l we can now derive the 
definition of the binding energy of an exciton in the presence of a localization potential: 

Eb{X, D) = Eb{X, 0) + <5^coui(X, D) + SEUX, D) + SEi^X, D), (14) 

where we define the change of the localization energy due to the excitonic bound state 

6EiUX,D) = 6EUX,D) + 6EtiX,D), 
^Et':\x,D) = Et':\D)-Et':\x,D) = j dr^^n)[9\r^iH))-9\r^iH).X)]Vi,"\ (15) 

Similarly, the change of the kinetic and interaction energies of an electron-hole pair forming an 
exciton is given by 

5EUX.D) = EUD) + EUD)-EUX,D) 

5EcoAX,D) = -ElUX,D). (16) 
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FIG. 5: Localization energy of an electron (a) and hole (b) in the different excitonic complexes as a function 
of diameter of the localization island for a one monolayer fluctuation and QW width of L = 60A. Dotted 
lines are the corresponding localization energies of a single electron and hole in the same locahzation 
potential. 

These expressions can be generalized directly to the case of the trions and biexciton. For example, 
for the positive trion the change in the Coulomb interaction is expressed as follows: 

5Ecoui{X\ D) = {Et.^,{X)} - {2E^cUX^) + ^coui(^+)}- (17) 

In Eqs. (fT6t - ([T7t the Coulomb energy, EcouiiD), is estimated as an average of the effective 
potential, r^g^ (Fig. H), over the pair distribution functions calculated for each type of interparticle 
interaction. For example, for the electron-hole interaction in the exciton we have: 

^coui(^) = jdr V:y{r) g,^{r, X). (18) 

The kinetic energy of the localized single electron (hole), E^l^\ the localized exciton, E]an{X), the 
trion, E]^a{X'^), and the biexciton, E]^a{X2), were computed as the difference between the total 
energy and the full potential energy which includes both Coulomb interaction and localization 
energy, E]^^ = E — {Ecovi + Eioc)- This result can be compared with a more strict thermodynamic 
estimator of the kinetic energy as the mass derivative of the partition function, E^^ ~ 
found that both expressions give very similar results. 
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In Fig.|5]we show the contribution of the well width fluctuation, which we will call the local- 
ization energies, to the single (non-interacting) particles, -Ef^^^ and interacting electrons (holes) in 
various excitonic complexes: K?^(^^)> K?^(^2)- This energy was estimated as an 

average of the localization potential, vj^'^'^^ (Eq. over the radial distribution functions calcu- 
lated separately for electrons and holes, see e.g. Figs.|7](c,d). For large D, E^^^^ should approach 
the depth of the in-plane localization potential, i.e. \V^'"^\ = 3.43meV and \Vl°''\ = 1.28meV. As 
one can see the interparticle interaction significantly increases the localization. Both, in the left 
and right panels of Fig. 13 the localization energy of a single particle (electron or hole) shown by 
the dotted line is much less (in absolute value) compared to that of the X, X^, X2 and saturates 
only for D > 1200 A. The explanation is that the attraction between electrons and holes already 
leads to significant spatial localization of the particles compared to the free particle thermal wave- 
length and, as a consequence, the effective localization potential felt by each particle (the potential 
is smoothened over the particle's wave functions) is deeper. 

An interesting point that we notice from Fig.|5]is that, on average, the electrons and holes in the 
biexciton are more localized than in all other bound states. The only exception is the localization 
energy of the hole in the range 200 A < D < 300 A when a single hole in the state tends to 
be more localized, see Fig.|5tb). Among all considered bound states only the biexciton appears to 
be strongly localized for islands with a diameter around D ^ 100 A (other excitonic states, in our 
simulations at temperature T ^ 1.5 K, have a much higher probability to become delocalized due 
to thermal fluctuations). This is confirmed by the gain in the binding energy shown in Fig.|3tb), 
where at the point D = 100 A only the biexciton shows a strong increase by 40%, and in Fig.|5]at 
D = 100 A only the biexciton shows nonnegligible values for the localization energy, Ei^^ = 0.54 
meV and Ej^^^ = 0.25 meV, for the electron and the hole respectively. 

Now, using the values of the localization energies of the electrons and the holes shown in 
Fig- El one can estimate the total localization energies of the different excitonic complexes. E.g. 
the localization energies of the negative trion and the biexciton can be estimated as follows: 

EiUX-) = 2EUX-) + E^X-), E1UX2) = 2EUX2) + 2EtiX2). (19) 

These energies include many-body correlation effects in combination with the specific radial dis- 
tribution functions which are different for each bound state, see e.g. Figs.|7](c,d). 

We now analyze the binding energies which are modified from their ideal expression, Eq. 
due to the three localization corrections, cf. Eqs. (I14I15I) . namely the contributions due to 
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FIG. 6: (a) Localization energy gain, 6Eioc (Eq. (flSl ). and energy gain, 6E, of the exciton and biexciton vs. 
diameter of the locaUzation island, (b) Comparison of the same energies for positive and negative trions. 
The same parameters are used as in Fig. 5. 

changes of the Coulomb interaction, kinetic and localization energy. Each part can be calcu- 
lated separately and the results are presented in Figs. |6ta,b). Solid lines in this figure show 
changes of the localization energy, 5Eioc{D), of the exction, trions, biexciton as a function of 
the diameter of the island, D. Dotted lines show combined changes in the Coulomb interac- 
tion and kinetic energy minus the binding energy in the same QW without the interface defect, 
5E{D) = SEcouiiD) + 5E^^{D) - ^b(O). Consequently, 

Eb{D) = Eb(0) + 6EUD) + SE{D). (20) 

From Fig. I^a) we note that the difference in the localization energies, shown by the solid 
curves in Figs. |6ta,b), is maximal for the exciton. This is easy to understand because for the 
exciton (see Eq. (fTSl) ) we subtract from the energies of the unbound electron and hole (which are 
less localized) the localization energy of a more localized bound electron-hole pair in the exciton 
state (see Fig.|5l). According to Eq. (l20b this gives positive contributions to the binding energy, e.g. 
at the point D = 300 A it is about 2 meV. For the biexciton and the negative trion the difference of 
localization energies reaches a maximum value of about 1 meV. For the biexciton the maximum 
is reached around the defect diameter D ^ 100 A, for X~ around D ^ 300 A, and for X+ at 
D ^ 200 A. 

The dotted lines in Figs.|6ta,b) are the contributions due to the defect induced changes of the 
Coulomb interaction and kinetic energy. In this way the total binding energy (l20b is the sum of the 
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FIG. 7: (a,b) Radial distribution function of electrons (a) and holes (b) for exciton, positive and negative 
trions, and the single electron (hole) as a function of distance, r^(^fj^y from the center of the localization 
potential. The same parameters are used as in Fig. 5. 

respective solid and dotted curves in Fig.|51 By comparing each pair of curves for X, and X2, 
one can note a general feature valid for all bound states: at D < 300 A the main contribution to 
the binding energy comes from changes of the localization energy, Eq. (fTSt . but when the defect 
diameter exceeds D ^ 600 A the main effect is due to changes of the kinetic energy and the 
interparticle interaction. 

It is interesting that for the exciton and biexciton the quantity AE becomes negative in the range 
of defect diameters D < 300 A which leads to a reduction of the binding energy. The reason is 
that 5E is composed of 5-Ekin which is negative and the positive term ^i^coui- ^-E'kin is negative 
because the kinetic energies of bound particles are larger than those of unbound particles. This 
is readily understood because the wave functions of bound particles are more localized (spatially 
less extended) and thus have a larger curvature, which increases the kinetic energy. It follows also 
from the virial theorem that the increase of potential (or interaction) energy leads to an increase of 
the kinetic energy. 

The comparison of SEioc and 5E between the positive and negative trions in Fig. |6tb) shows 
that the negative trion is more affected by the localization for D > 100 A. This has a direct relation 
with Fig. 13 where the X^ binding energy exceeds the one of the X+ for all D > 150 A. 

The positive trion X"^ is a much heavier composite particle with two holes and one electron 
compared to the X^, where there are two electrons and one hole. As a consequence the X+ is 
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less mobile and, thus, is less affected by the lateral confinement. This fact can be demonstrated 
by Figs. |7ta,b) where we show the radial distribution functions g^{re(h)) of an exciton, nega- 
tive/positive trions and a single electron/hole in a QW of width L = 60 A for a defect diameter 
D = 300 A. The panels (a),(b) show the radial density of electrons and holes, respectively. Notice 
that at the center of the localization confinement (re(/i) = 0) the radial density of electrons in the 
state is increased by more than a factor of 3 as compared to the radial density of a single 
electron in the same localization potential (curve indicated by "e"). At the same time, compari- 
son of the densities for the holes in the X+ state and the density of the single hole shows only 
an increase by a factor of 1.8. These results are in agreement with our previous statement that 
the localization has a larger effect on electrons than on holes and, consequently, on the composite 
particles with a larger number of electrons (in our case and X2). From Figs. |7|;a,b) we can 
also see that the central radial density of electrons is highest in the X+ state, whereas the highest 
central radial density of holes is observed in the X^ . In both panels (a) and (b) the electron and 
hole radial densities near the localization potential center for the exciton state lies in between the 
corresponding values for the X+ and X^. This fact is easily understood from the symmetry of the 
spatial configuration of particles in the X, X+ and X^ states in the cylindrical defect potential. 
For trions, the single hole in the X^ (the single electron in the X+) will most probably occupy the 
central position in between the two electrons (two holes) to minimize the correlation energy and 
thus the total energy. 

C. Dependence on the number of monolayer fluctuations 

Now we allow for well width fluctuations larger than 1 ML and analyze the dependence of the 
binding energies and interparticle distances on the number of monolayers. We fix the quantum 
well width to L = 40 A and the diameter of the defect to D = 400 A. In Fig. [S] we plot the binding 
energy of the different excitonic complexes as a function of the number of monolayers forming 
the defect (the curves are guides to the eye since is a discrete index). Notice that increasing 
leads to a monotonic increase in the binding energy for all exciton complexes which saturates for 

^ 4. Notice that the increase of the binding energy of the biexciton and of the negative trion is 
almost parallel. This is a clear confirmation of our earlier conclusion that the lateral confinement 
of the electron by the defect has a more pronounced effect on the binding energy than the hole 
confinement. Again we observe that the binding energy of the X~ exceeds that of the X+ in the 
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FIG. 8: The same as Fig. 13 but now as a function of the number of monolayers, N, for a fixed QW width 
of L = 40 A, temperature T = 1.5 K and diameter of the localization potential, D = 400 A. 
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FIG. 9: The same as Fig.|4j but now as a function of the number of monolayers, A'^. Same parameters are 
used as in Fig.[Sl 

presence of a localization potential (see the discussion of Fig. |3l). Fig. [H shows that this trend 
persists in the case of increasing defect depth. 

Fig.|9ldisplays the dependence of the mean e-e, h-h and e-h distances in the different excitonic 
states on the number of monolayers. All distances (i.e. the spatial extension of all bound states) 
decrease monotonically with N and saturate around = 4. 
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V. COMPARISON WITH EXPERIMENT 



Before making a comparison between our theoretical results and the experimental data it is 
necessary to point out a problem that may arise in such a comparison. In particular, one has to 
be aware of the fact that for finite temperatures there exists an experimental uncertainty about the 
state in which particles remain after recombination. Namely, particles produced after recombi- 
nation will typically have a finite kinetic energy which may allow them to leave the localization 
potential. In this case, the energy of the emitted photon will be reduced by an amount needed to 
overcome the height of the localization potential. As a consequence, the photoluminescence lines 
of the exciton complexes exhibit an additional broadening due to the finite kinetic energy of the 
remaining particles, thus making the determination of the binding energy more complicated. 

Secondly, in experiments there may be two types of measurements of excitons and trions. One, 
when the excitonic complexes are probed in a single well. These measurements may favor ob- 
servation of the most strongly localized excitons and trions (as it was found in Ref. |29]), while 
measurements from a QW ensemble favor higher-lying states nearer to the continuum. In the lat- 
ter case the localization effects does not strongly influence the observed ensemble QW spectra. 
As the theory shows, this would affect the binding energy, as the localization effects are of great 
importance in narrow QW's. 

A. Exciton binding energy 

We first compare in Figs. fTOl' a.b') the theoretical and experimental binding energies of the exci- 
ton as function of the QW width. Solid and dashed curves in these figures show the binding energy 
in the QW with defect and in the ideal QW without interface roughness, respectively. The local- 
ization potential is considered as due to a well width fluctuation of one monolayer over a circular 



area of diameter D = 400 A in accordance with the experimental findings of Ref. U^. As we can 
see from Fig.|3]for the QW width around 60 A this gives an upper bound to the localization effect 
on the binding energies. In order to be consistent we used the same localization potential when 
calculating the binding energies of trions and the biexciton. 

When using the isotropic approximation for the hole mass (mg = 0.067 mo, rri^ = = 
0.34 mo), we notice from Fig. flOl' a') that the theoretical results systematically overestimates the 
binding energy. For example, for L > 100 A the theory gives binding energies which are approx- 
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FIG. 10: Exciton binding energies for isotropic (a) and anisotropic (b) hole mass vs. quantum well width. 
Solid square symbols are experimental data of Ref. 1^ and circular symbols with a plus sign are obtained 



from the stochastic variational calculations 



141 . Full (dashed) lines with small symbols are PIMC 



results when the localization is included (not included). 



imately 10% larger than the experimental results. In order to see the effect of the anisotropy of 
the hole band we performed the calculation using an anisotropic hole mass (m^^ = 0.112mo in 
the QW plane and = 0.377mo in the QW growth direction [26]) which brings the theoreti- 
cal points significantly closer to the experiment, see Fig. fTOl' b). For the QW widths L < 100 A 
good agreement is found when the localization effects are included in our model. For example, in 
Ref. [3] for a QW width L = 34 A a value 12.4 meV has been measured for a localized exciton. 
This agrees quite well with our theoretical prediction of about 12.3 meV for a comparable QW 
structure (L = 40 A). 

However, there is still a slight discrepancy between theory and experiment for QW widths in 
the range 100 < L < 200 A. In this range, in our theoretical model, the localization effects are 
negligible for excitons. For L > 150 A the two curves calculated with and without localization 
practically coincide. In this case, for the wide QW's (L > 150 A) the theory appears to agree quite 
well with the experimental binding energies. This allows us to conclude that for such quantum well 
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structures the localization does not play significant role for excitons and they are not trapped by 
the interface defects. 

In Figs. [TUt a.b) we compare our results with those obtained with the stochastic variational 
approach jl3[|l4l . both for the isotropic and anisotropic hole masses which gives additional credit 
to the accuracy of our numerical approach. Note that the hole mass and the dielectric constant 



racy or ou 



used in Ref. II13L 11411 (me = 0.067mo, = 0.099 mo, e = 12.1) were slightly different from 
ours. This leads to minor discrepancies for the binding energies for QW widths smaller than 100 
A because the electron and hole densities in the growth direction are very sensitive to the QW 
confinement and to the chosen values of electron/ hole masses in narrow QW's (see Fig.|2ll. 

The PIMC results (for the anisotropic hole mass and the 1 ML interface defect of diameter 



2lh . where the exciton binding energy 



D = 400 A) can be compared with the those of Refs. 112(1 
to the interface defect, E^{X) has been calculated. This quantity is defined as the difference 
of the total exciton energy with and without the defect potential. In particular, the variational 
calculations of Ref. [21] for two types of exciton trial wavefunctions give, for a QW of width 
L = 35 A and D = 400 A, E^{X) = 7.4 meV and 10.4 meV, respectively (the isotropic electron 
mass is ml = m^^ = 0.0782mo, the heavy hole masses are the same as in the present work). The 
PIMC calculations for the QW width L = 40 A give a similar value, E^{X) = 8.16 meV. In the 
variational calculations of Ref. a value E^{X) = 3.4 meV has been reported for a QW with 
L = 70 A. The present calculations give E^{X) = 3.72 meV and E^{X) = 0.84 meV, for a QW 
width L = 60 A and L = 100 A respectively. This trend definitely shows that the excitons become 
less localized with increasing QW width. 

In conclusion, the comparison of Figs. flOl' a') and (b) shows that taking into account the 
anisotropy of the hole mass leads to a decrease of the exciton binding energy compared to the 
isotropic case. For example, in a 50 A wide QW this amounts to about 2 meV, and in a 250 A wide 
QW it is about 1 meV. 



B. Binding energy of positive and negative trions 



In Figs. fTTT' a.b') we present our results for the binding energy of the trions. We compare our 
theoretical results with the available experimental data for negative 1 1| and positive [Q] trions, and 
variational calculations of Ref. ll22ll . First, we can note that the agreement with the experiments 
is quite good for QW widths L > 150 A. Specifically, the experimental points for the X~^ [j2|] are 



23 




' — I — I — I — ' — r 

-0 — isotropic, X 

isotropic, X* 

* X " Shields ef a/. 
■ X Kaurefa/. 
A X Finkelstein et al. 

• X Yanefa/. 
X * Osborne et al. 

o X* Blakemore et al'. 
< X* Phillips ef a/. - 



100 200 300 400 
L(A) 







4 
3 
3 
2 
2 
1 
1 



100 200 300 400 
L(A) 



1 ' 1 
- 


— — anisotropic, X . 
— o — ani^otrnnir 
V X * Dacal ef. al. 
(anisotropic) 






. " - V 




: b) 
1 . 1 







5 


5 "o 

- 2.0 ^r- 

+ 

5^ 

LU 



5 




FIG. 11: Trion binding energies for isotropic (a) and anisotropic (b) hole mass vs. quantum well width. 

\^\^ and theoretical calculations of Ref. I^l . 



Symbols are experimental data of Refs. ll|,|2|] and theoretical calculations of Ref. 
with symbols are PIMC results when localization is included (not included) 



|. Full (dashed) lines 



close to our theoretical result, see Fig. fTlT a). Unfortunately, for narrow QW's (L < 150 A) when 
localization effects become important, there are currently no available experimental data. This 
would be of high interest as the two points for the X reported by R. Kaur et al. and Z.C. Yan et 
al. Ill show a more rapid increase of the binding energy with the QW width than the one predicted 
by theory when the localization is not taken into account. On the contrary, calculations with the 
QW width fluctuations included agree well with these data. 

Further, in Fig. fTlT a) we observe a crossing of the two binding energy curves (X+ and X^) 
with the localization included near the point L = 135 A, and for narrow wells the binding energy 
of the negative trion becomes larger. For example, in a 40 A quantum well Eb{X~) ^ 3.4 meV 
while for the positive trion Eb{X^) ^ 2.6 meV. Preliminary experimental investigations of the 
X~ and X^ in the presence of localization ll29ll seem to confirm these findings and give for a 50 
A wide QW Eb{X-) = 3.27 meV and Eb{X^) = 2.35 meV. It is interesting to notice that in an 
ideal QW the situation is the opposite and the binding energy of the X+ is always larger than the 
one of the X . As was already discussed in section HVBl the reason is that the localization has a 
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stronger influence on the electrons and, consequently, on the negative trion. 

When the anisotropy of the hole mass is included in our calculations we found that there is no 
crossing between the binding energy of the X+ and the . A comparison with the variational 
calculations of L. Dacal et al. |22ll (see Fig. [HJb)) also done with the anisotropic hole mass, 
but slightly different parameters VeiVh) = 224.5 meV (149.6 meV) and e = 13.2 (compared to 
ours Ve{Vh) = 216 meV (163 meV) and e = 12.58), shows some discrepancy with the present 
calculations for X+ binding energy in the narrow QW's. This can be due to the following reason. 
Variational calculations strongly depend on a form of used trial wave functions. In particular, the 
Eq. (3) in Ref. Q becomes less accurate in the narrow quantum wells and hence requires a larger 
set of variational parameters. This could lead to a better agreement with the PIMC results. 

In Ref. ll23h the effect of localization (at the interface defect with a cylindrical symmetry and a 
Gaussian shape) on the trion binding energy was considered for the anisotropic hole mass with the 
same parameters as in Ref. 11221] . It was found that, with the defect present, the X" binding energy 
is increased from 0.4 to 0.6 meV in the 150 A wide QW. However, these values are much lower 
than the PIMC results which show an increase from 0.90 to 1.02 meV. The results of Ref. ll23l] 
are even lower than the value Eb{X^) = 0.75 meV reported in Ref. for the same QW 
width but without the localization effect included. In Ref. ll23h the number of basic variational 



trial wavefunctions was reduced compared to Ref. 112211 and, as the above comparison shows, this 
appears to be not sufficient for a quantitative description of the trion. 

Other theoretical calculations done with an isotropic hole mass and in the absence of a local- 
ization potential shown in Fig. fTTV a) agree quite well with our data for the ideal QW case. For 
example, both in Ref. 1 2Q ] and Ref. Il6ll it was found that the X+ binding energy is larger by about 
20% than the binding energy. This is in agreement with most experimental results which show 
that the X^ has a binding energy which is larger than or close to the one of the X~ [1]. How- 
ever, the theoretical results of Ref. Ini for a 300 A wide QW showed that Eb{X^) is lower than 
Eb{X^) which is opposite to the results of Refs. Il^ Isoi . but the latter is in agreement with 
our results for the case of the anisotropic hole mass (see Fig. 11(b)). From the other hand, the 
anisotropic calculations show not satisfactory agreement between theory and experiment, in par- 
ticular for L > 150 A where theoretical curves are about 0.5 meV below the experimental points. 
Here the localization effects (in the framework of our model) are not important and can not be the 
reason for this disagreement. 
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C. Binding energy of biexciton and the Haynes factor 



Now we consider the biexciton binding energy. Solid (dashed) curves in Figs. [T^ a.b) show 
our theoretical results in the presence (without) of well width fluctuations which we compare with 
available experiments lOiyiyiyiyiSQ] (symbols in figure). In the case of an isotropic hole mass, 
we found that with localization taken into account (full line in Fig. [12 a)) the theoretical curve 
passes through the data of Refs. ijs.Q], but the experimental data of Refs. 10,1^ agrees with 
the theory for an ideal QW (see the dashed curve in Fig.[T3a)). The inclusion of the localization 
effects brings the theoretical curve slightly above the experiment, predicting that in e.g. a narrow 
40 A wide QW the binding energy of the biexciton is about(4.0 — 4.1) meV. It is interesting that 
practically the same experimental result is reported in Ref. llOll . where the value 4. 1 meV is found 
for the 40 A QW. It should be stressed however, that our results for the binding energy, in the 
presence of a one monolayer well width fluctuation, are close to an upper limit, since the used 
defect diameter (D = 400 A) was such that it gave practically the maximal gain in the biexciton 
binding energy as due to localization (see Fig.|3l). But this value of the defect diameter, D = 400 A 
was found to be a very good estimate of the characteristic defect size in GaAs quantum wells liioi . 
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The use of an anisotropic hole mass (see Fig. [T2r b)) leads to a reduction of the biexciton binding 
energy by almost (0.4 — 0.8) meV. As in the case of trions, for wide QW's with L > 150 A the 
anisotropic approximation gives a less satisfactory agreement with the experimental points. From 
the other hand, for the QW widths L < 150 A we found marvellous agreement with the experiment 
ofRef. 

when localization is included in our calculations (see Fig. fT^ b). solid curve). 

For an anisotropic hole mass and a defect depth of 1 ML, we can compare our results with the 
variational calculations of Ref. l2li . In these calculations, however, a repulsion between particles 
of the same charge was not taken into account, which makes the comparison only qualitative. In 
Ref. [Illl the biexciton binding energy of about 1.6 meV and 1.3 meV have been reported for the 
30 A and 50 A wide QW's, respectively, and a defect diameter 200 A. These values of the binding 
energies are less than our result, 1.9 meV, for a non-localized biexciton in a 50 wide A QW. 
With the interface defect (D = 400 A), the PIMC calculations show an increase of the binding 
energy up to 3.6 meV, which is close to the experimental data. For example, in Ref. [^, a value 
Eb{X2) = 4.2 meV was attributed to a localized biexciton in a 34 A wide QW. 

In conclusion, the present comparison of the theory and the experiment allows us to conclude 
that the anisotropic approximation for the hole mass gives better agreement for narrow QW's 
(with L i 150 A). On the other hand, for wide QW's (L > 150 A) our model calculations show 
that the use of isotropic approximation show better agreement. We expect that the accuracy of 
the calculations (and of the anisotropic approximation, in particular) can be further improved by 
taking into account the mismatch of dielectric constants and particle masses in the well and barrier 
materials. This should lead to a better agreement with the experiment. 

In comparing the theory and the experiment one should keep in mind that different experimen- 
tal results have been obtained from different quantum wells which have not been grown under 
the same conditions and, consequently, their well width fluctuations may also be different. Nev- 
ertheless, in overall, the present calculations show that even a simple model of localization can 
satisfactory explain the experimental data on the binding energies. We expect that our theory can 
make reliable predictions if detailed information on the quantum well fluctuations in the samples 
would be available. 

Next, we consider the so-called Haynes factor, which is the ratio between the biexciton and the 
exciton binding energies, u = Eb{X2) / Eb {X) . We compare our results, see Fig. flBl a). with the 
experimental data of Ref. Iioi (solid squares), where the influence of localization on the binding 
energies of the exciton and biexciton was considered. Here all theoretical curves are given only 
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FIG. 13: The biexciton to exciton binding energy ratio (Haynes factor), (a) versus QW width and (b) vs. 

n 

normalized localization potential. Symbols: experimental data for GaAs/Alo.sGao.rAs QWs 1 10]. Curves 
with symbols: PIMC results. Solid (dashed) curve is the exciton binding energy with (without) localization 
included, (b) Eb{X2)/Eb{X) versus localization strength E\oc/ Eb{X2). Squares are the experimental 
data. Rhombus are PIMC results. 



for the isotropic hole mass (m^ 



xy 



mr 



0.34) and, as in previous figures, the solid (dashed) 



line is for the case (without) localization included. The ratio of the exciton and biexciton energies 
calculated using an anisotropic hole mass is very similar and is therefore not shown. 

The lines with small filled symbols are obtained by using the theoretical values of the exciton 
and biexciton binding energies (Figs. fTOl' a) and [I2a)), which are denoted as E^^^{X2)/ E^^'^{X) 
and Eb{X2) / Eb{X). In addition, since there is a discrepancy between the experimental and the 
theoretical results in Fig. fTOl a) for the exciton binding energy which affects the Haynes factor, 
we calculate also the Haynes factor as the ratio of the theoretical biexciton binding energy to 
the exciton energy from the experiment lioi (shown in Fig. fTOl a) by the solid squares). The 
corresponding results are denoted in Fig.[l3l;a) as E^^''{X2) / Esy,p{X) and Eb{X2)/ Es^p{X). 

First, we note that the Haynes factor is practically independent of the QW width for L > 
150 A in agreement with experiment. However, the value of the constant is different in various 
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cases. When localization is included (solid line) our calculations show a systematic increase of 
the Haynes factor v = Eb{X2) / Eb{X) by up to 17%. For the exciton and biexciton localized 
on the interface defect the Haynes factor is z/ ^ 0.175, while for the ideal QW, v ~ 0.15. For 
well widths L > 50 A, all theoretical and experimental results exhibit a monotonic decrease of the 
Haynes factor with increasing L. One can observe that the experimental results are located mainly 
between the two theoretical curves corresponding to the localized and not localized biexciton. 
Most of the points lie on the dashed curve with the open circles, Eb{X2) / E^t^^^X), suggesting 
that the experimentally measured binding energies correspond to non-localized biexcitons, but 
some of the data are substantially above the dashed curve and agree better with the assumption 
of predominantly localized biexcitons. This is the case for well width L < 80 A where a strong 
increase of the Haynes factor is found also from the theory. For example, our theory gives a 
maximum value u = 0.24 without localization and u = 0.29 for localized particles. 

Notice that the agreement between the experimental and theoretical results with localization is 
also confirmed by Fig. fTSl b). where the Haynes factor is plotted against the normalized localization 
potential. We compute the localization energy Eioc as the difference between the energy of the 
localized and non-localized excitonic complex (in Ref. |Q] the localization energy was defined 
from the full width at half maximum of the heavy-hole exciton absorption line). As follows from 
our theory, below 150 A well thickness when the localization energy becomes of the order of 
the biexciton binding energy an enhancement of the Haynes factor is observed. This indicates that 
localization has a crucial effect on the Haynes factor and it must be taken into account for a correct 
interpretation of the experimental results. 



VI. CONCLUSION 



In summary, using a path integral Monte Carlo approach, we made a detailed analysis of dif- 
ferent excitonic complexes in GaAs/AlGaAs quantum wells. We calculated and analyzed the 
exciton, trion and biexciton binding energies, pair distribution functions and mean interparticle 
distances in the excitonic complexes in a wide range of QW widths. Our method is based on first 
principles and does not invoke expansions in eigenf unctions. The approach is general, flexible 
and - what is also important - is not limited to certain specific symmetries of the particle wave 
function. It allows for a simultaneous account of QW confinement, localization and valence band 
anisotropy, and thus can give an accurate theoretical treatment of many experimental systems. The 
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only assumption - the adiabatic approximation (Eq. Q) - appears to be justified for the present 
application. Simple estimates show that its accuracy may be reduced for wide QW's with L > 250 

A. 

Extending our previous analysis I24I, we concentrated here on the influence of disorder, i.e. 
of the effect of QW width fluctuations on the binding energies. The observed increase of the 
binding energies in the presence of disorder is in good quantitative agreement with the available 
experimental data. Furthermore, we analyzed the influence of valence band anisotropy (hole mass) 
and found that, in some cases (in particular for the exciton and biexciton binding energy) this effect 
is important in order to achieve agreement with the experimental data. 

We also analyzed the case of deep interface defects corresponding to several monolayers depth 
and found that this can give an additional 20 — 30% increase of the binding energy compared to 
the single monolayer case. This increase is even more pronounced for the negative trion and the 
biexciton. 

The present analysis is the first one in which exciton, trions and biexciton are treated on an 
equal footing and in which the same size and shape of the QW width fluctuation is invoked for 
different QW widths. We assumed a IML QW width fluctuation over a circular area of diameter 
D = 400A. No other fitting parameters were introduced. This lead to an overall good agreement of 
the well width dependence of the exciton, trions and biexciton binding energies. It is expected that 
a better fit with experiment is possible if, e.g. we allow for a non circular shape of the well width 
fluctuation where the anisotropic localization potential will be related to specific crystallographic 
directions. 

Our results have two important implications which can be useful in the interpretation of exper- 
imental data. First, by comparing the measured binding energy with our numerical calculations 
for different defect sizes allows one to characterize certain experimental parameters, such as the 
magnitude of the disorder in a given sample. Secondly, one can verify or predict whether or not 
the observed excitonic states are localized or delocalized in a given experimental set up. 
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